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ABSTRACT 

We study constraints from the latest CMB, large scale structure (2dF, Abell/ACO, 
PSCz) and SNla data on dark energy models with a sharp transition in their equation 
of state, w(z). Such a transition is motivated by models like vacuum metamorphosis 
where non-perturbative quantum effects are important at late times. We allow the 
transition to occur at a specific redshift, Zt, to a final negative pressure — 1 < Wf < 
— 1/3. We find that the CMB and supernovae data, in particular, prefer a late-time 
transition due to the associated delay in cosmic acceleration. The best fits (±l<r errors) 



to all the data are Zt = 2.0 



+2.2 
-0.76' "Q 



0.73 



+0.02 
-0.04 



and Wf = — 1+° • 



1 INTRODUCTION 



The idea that the universe is currently accelerating comes 
from a number of high-quality, but indirect, experiments. 
The luminosity distance estimated from type la supernovae 
(Riess et al. 1998) favours recent acceleration while the re- 
cent Cosmic Microwave Background (CMB) data (Netter- 
field et al. 2002, Lee et al. 2001, Halverson et al. 2002, 
Padin et al. 2001) suggest the universe has almost zero 
spatial curvature. This, combined with clustering estimates 
giving Q m ~ 0.3 (Hamilton and Tegmark 2002), provides 
compelling evidence for a dominant, unclustered, univer- 
sal element; a conclusion supported by the height of the 
first, and position of the second, acoustic peak in the CMB 
(Kamionkowski & Buchalter 2000). 

Nevertheless, our understanding of the true nature of 
such an unclustered "ether" is arguably even worse than it 
is for the dark matter responsible for galaxy clustering. The 
oldest idea is that it is a cosmological constant, A^0. This 
requires generating a tiny scale A ~ (10 _3 eV) 4 and there 
have been many recent ideas to achieve this (see e.g. Dienes 
2001, Deffayet et al. 2001, Verlinde & Verlinde 2000, Ver- 
linde 2000), mostly motivated by the revolution associated 
with higher dimension brane world models. 

Undoubtedly the best-studied explanation, however, is 
quintessence (Ferreira & Joyce 1998, Caldwell et al 1998); a 
very light scalar field Q, whose effective potential V(Q) leads 
to acceleration in the late universe. However, quintessence 
suffers from extreme fine-tuning since not only must one set 
the cosmological constant to zero but one must arrange for 
the quintessence field to dominate at late times only. This co- 
incidence problem typically requires severe fine-tuning in the 
potential, although this can be alleviated in models where 
the quintessence field couples to dark matter (Holden & 
Wands 2000, Amendola 2001, Tocchini-Valentini & Amen- 
dola, 2001). 

Such couplings to standard model fields imply ex- 



tra fine-tunings however since the quintessence field is ex- 
traordinarily light and typically has Planck-scale vacuum 
expectation value today Q ~ M v i. If one believes that 
the quintessence field comes from supergravity, then non- 
renormalisable couplings to standard matter fields should 
appear automatically, even if all the renormalisable cou- 
plings to standard matter are fine-tuned to vanish. 

These non- 
renormalisable couplings, such as QF^F^" /M p i, where F^ v 
is the Maxwell tensor, cause variations of the fundamental 
constants of nature, such as the fine-structure constant (Car- 
roll 1998), and the corresponding coefficients must be made 
small to avoid conflicting with evidence at low redshift. Sim- 
ilar constraints come from neutrino data (Horvat 2000) and 
make the quintessence scenario appear extremely fine-tuned 
without some more fundamental theoretical basis. 

Quite another possibility is that quantum effects have 
become important at low redshifts and have stimulated 
the universe to begin accelerating. Examples are vacuum 
metamorphosis, put forward recently by Parker & Raval 
(1999,2000: hereafter PR) and the work of Sahni and Habib 
(1998). 

In particular, PR consider a scalar field of mass m in 
a flat, FLRW background, and non-perturbatively compute 
the effective action in terms of the Ricci scalar, R. They show 
that the trace of the Einstein equations receives quantum 
corrections some of which are proportional to 



hGm 4 R 



m? + (£ - 1/6)7? 



[1 + 0(R)]. 



which diverges when R 



(1) 



7(£ - f/6), signalling sig- 



nificant quantum contributions to the equation of state of 
the scalar field. Here £ 7^ 1/6 is the non-minimal coupling 
constant. At early times the equation of state is dust on aver- 
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age [f] and then makes a transition from dust to cosmological 
constant plus radiation (PR 1999). 

To explain the supernova Type la (SNla) data the 
scalar field is forced to be extremely light, m 2 /(£ — 1/6) ~ 
10 -33 eV. Vacuum metamorphosis therefore suffers from the 
same fine-tuning problems as quintessence - why are there 
no dimension 5 operators leading to unacceptable variation 
of fundamental constants? 

The idea of a sudden phase transition is very attrac- 
tive however and is more general than just the example of 
vacuum metamorphosis. In fact the idea of late-time phase 
transitions is rather old, dating back at least as far as 1989 
(Hill et al 1989, Press el al 1990). 

We therefore choose a phenomenological model which 
captures the basic features of a phase transition in the equa- 
tion of state, but which is not strictly linked to any specific 
model. We then ask whether current CMB and large scale 
structure (LSS) data rule out such a transition, or indeed, 
favour it over the now standard ACDM model. 



2 THE PHENOMENOLOGICAL MODEL 

In addition to baryons, neutrinos and cold dark matter our 
model is characterized by a scalar field Q with a redshift 
dependent equation of state pq = w(z)pq. We choose w(z) 
to have the following form 



w(z) = Wo + 



(Wf — w ) 
l + exp(^ 



(2) 



In this paper, we shall restrict ourselves to the case where the 
initial equation of state is pressure-free matter, wo = 0. In 
terms of the scalar field dynamics, we will consider a class of 
models characterized by three free parameters: (1) the final 
value of the equation of state, Wf, (2) the redshift z t of the 
transition, and (3) the energy density Hq of the scalar field 

in units of the critical energy density p cr it', Qq = Pq/Po-u f] 
We also assume that any coupling of the scalar field with 
other fields are negligible. In this case the energy density pq 
is determined from energy conservation 



Pq = ~3Hp Q (l + w(z)) 



(3) 



which can be explicitly integrated since w(z) and Q,q today 
are given. Using Eq. (bl) for w(z) and specifying initial condi- 
tions for the scalar field one obtains the scalar field potential 
V(Q) and its derivatives along the "background" trajectory 
Q(t). In particular, V'(Q), is 



i_i ( i +z) dl °gK*)) 

3 dz 



(4) 



t A massive scalar field acts like dust on averaging over many 
oscillations of the field. Since the scalar field here is so light this 
may be a bad approximation since the period of oscillation is so 
long. 

t Notice that the width of the transition is controlled by A. Dou- 
ble precision limits require that zt/A < 100 while the constraint 
that w = wjahz = § implies that zt/A > 10. To satisfy both 
of these constraints simultaneously for < zt < 2000, we have 
chosen a zt-dependent A defined by zy/A = 30. 



and V"{Q) can be easily obtained from Eq. (Q). V"(Q) is 
required to solve the evolution equation for the scalar field 
fluctuations, <5Q(x, t). In the synchronous gauge, the Fourier 
modes SQ^ (t) obey the equation 



h + 3HSQ k + [ !L + V" 



-h k Q 



(5) 



where h is the trace of the spatial metric perturbation (Ma 
and Bertschinger 1995). We choose adiabatic initial condi- 
tions 5Qh = SQk = for the scalar field. 



3 THE DATA AND ANALYSIS PIPELINE 
3.1 The cosmic parameters 

Due to computational restrictions we fixed cosmic parame- 
ters not directly linked with the scalar field Q. We therefore 
performed a likelihood analysis in the neighbourhood of the 
best fit standard model. Hence we have no assurance that a 
better global minimum for does not exist. However, our 
main goal is to test whether the current data favour a phase 
transition over the standard ACDM model, and for this our 
analysis is sufficient. We chose the following "standard" cos- 
mic parameters (Wang et al. 2001): 

• Ho = 65 km/s/Mpc 

• fib — 0.05, fitot = 1 (flat universe) 



• n. 



= 1, n t = 0, r = 0. 



Here r is the reionisation optical depth and n 3t t are the 
spectral indices for scalar and tensor perturbations respec- 
tively. We set Qa = and included only the effective 3.04 
massless standard model neutrinos. We did not include ten- 
sor perturbations and as we varied SIq today we specified 
the cold dark matter density to ensure overall flatness of the 
universe, viz. SIcdm = 1 — fit — Qq. We emphasize here that 
fixing fib and Ho can produce artificially narrow likelihood 
curves especially for Qq and that our results should not be 
interpreted as determining the true energy density of the 
scalar field to such a precision. 



3.2 Observational Data and Analysis 

We constrain the parameters of our phenomenological model 
by comparing its predictions with a number of observations. 

CMB : we use the decorrelated COBE DMR data 
(Tegmark & Hamilton 1997) to constrain the fluctuations 
on large scales and combine it with the recent data from the 
BOOMERanG (Netterfield et al. 2002), MAXIMA (Lee et 
al. , 2001) and DASI (Halverson et al. 2002) experiments. 
In total we used 49 data points ranging in £ £ [2, 1235]. We 
take into account the published calibration uncertainties for 
BOOMERanG, MAXIMA and DASI but not the pointing 
and beam uncertainties, since they are not public for all 
experiments. 

LSS : we use the matter power spectrum inferred from 
the 2dF 100k redshift survey (Tegmark, Hamilton & Xu 
2001), the IRAS PSCz 0.6 Jy survey (Hamilton & Tegmark, 
2000) and the Abell/ACO cluster survey (Miller et al. 2001). 
We limit the comparison to k < 0.2fa/Mpc to minimise po- 
tential non-linear contaminations. All together, we use 48 
points. We do not currently include the Ly-a analysis of 
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Croft et al (2000). Even though our best fitting models 
agree quite well with the shape of the recovered linear mat- 
ter power spectrum, the concerns raised by (Zaldarriaga et 
al. 2001) make it seem preferable to postpone the use of this 
data set. 

SNla : we use the redshift-binned supernova data from 
Riess et al. (2001), which includes the HZT (Riess et 
al. 1998) and SCP (Perlmutter et al. 1999) data. 

We follow the standard approach of computing the CV's 
and P(k) for each set of parameters over the 3 dimensional 
grid (flQ,ZT,Wf), using a modified version of CMBFAST 
(Seljak & Zaldarriaga 1996). We then evaluate the corre- 
sponding \ 2 a t each grid point. 

We find that the likelihood values for the CMB depend 
slightly on the likelihood functional used, eg. a simple \ 2 
computed using the CVs or an offset log- normal distribu- 
tion (Bond et al. 2000), but we are not certain if this is an 
intrinsic difference, or due to the slightly different data in 
the RADPACK package. However the overall changes are 
within la and hence not significant. 

The link between CMB and LSS is given by the respec- 
tive normalisations. The CMB data fixes the overall ampli- 
tude of the model quite precisely. As the connection between 
the matter power spectrum inferred from galaxy and cluster 
surveys and the actual distribution of dark matter is much 
less clear, we allowed a bias, b £ (1/5,5) for 2dF and PSCz 
and b £ (1/9,9) for Abell/ACO since clusters are expected 
to be more biased than galaxies. Here b is the factor between 
the perturbation amplitudes, and hence enters quadratically 
in the power spectra. 

In general we marginalise over parameters by integrat- 
ing the likelihood. We find that the results are consistent 
with those found from maximising the likelihood. This is 
expected for a nearly Gaussian likelihood, and it provides 
some reassurance that the \ 2 method is justified. 

4 THE PHYSICS OF METAMORPHOSIS 
4.1 The CMB 

To understand the imprint of metamorphosis on the Ci 's of 
the CMB requires two insights. First the contribution of the 
scalar field to the expansion rate of the universe is negligible 
for z > 3 if Wf < -0.4 since p Q ~ a~ 3( - 1+w f \ This implies 
that the dynamics of Q has little effect on the evolution of 
the metric perturbations (which respond to the total matter 
perturbation) and hence the CVs are almost insensitive to 
transitions with z t > 3. 

This is evident in Fig. (^). The figure also shows an ef- 
fect which at first sight is perhaps surprising: the CMB is 
extremely sensitive to Zt for z t < 3. This is clarified once we 
remember that the standard ACDM CMB has a large Inte- 
grated Sachs- Wolfe (ISW) contribution due to the decay of 
the gravitational potential $ during the epoch of accelera- 
tion, typically occurring at z < 2. 

Therefore, if we choose two models with wf < —0.5 but 
with zt = 0.5 and z t = 1.5 respectively, the initiation of 
the accelerated phase will vary by around 300%, and the 
decay of the gravitational potential starts at very different 
epochs. Hence the fluctuations in the CMB for these two 
models differ mainly on large angular scales which alters 
the acoustic-peak/SW-plateau ratio. 




1 k [h/Mpc] 

Figure 1. Left: The total-data best-fitting model (right curve) vs 
ACDM (left curve) which both have Uq = 0.73 and Wf = — 1 but 
zt = 1.5 for the best-fit. Right:The power spectrum for our total 
best-fit model with (Qq , zt,w f) = (0.73, 1.5,-1) compared with 
the LSS best-fit (0.7,6.5,-0.55). The LSS data shown are the 
linear transfer functions inferred from the 2df and PSCz galaxy 
surveys and the Abell/ACO cluster survey. We do not show the 
lyman-c* data. 
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Figure 2. Variation of the CVs with zt and Wf. Left: The Ci 
curves increase monotonically with decreasing Wf starting with 
-0.6 (bottom), -0.7, -0.8, -0.95 and ending with -1 (top). The 
primary change is in the normalisation of the spectrum through 
the change in the ISW effect. Right: The CVs for zt = 0.5 (bot- 
tom), 1.5,3,5 and 10. The ISW contribution to the COBE nor- 
malisation changes very rapidly for small zt which allows for de- 
layed acceleration. However, the CMB is insensitive to zt > 3 
since the scalar field is dynamically irrelevant at those redshifts. 

In addition, z t < 2 implies that the rapid change in 
w can cause its own effect on the gravitational potential 
since the energy density of Q is starting to dominate. This 
gives rise to an ISW effect purely due to the scalar field 
dynamics. We will see later that these two effects almost 
completely explain the behaviour of the likelihood curves in 
Fig.M) which are very sensitive to small Zt but exhibit a 
long tail for large Zt- 

The other parameter of interest is Wf. This has a sim- 
ple effect on the CMB. As w/ decreases towards to —1 the 
universe starts to accelerate earlier and is accelerating more 
violently today. This alters the ISW effect which is impor- 
tant on large scales and which contributes to the COBE 
normalisation. The effect of Wf is mainly then to amplify 
or supress the £ > 50 CVs by a more-or-less £- independent 
amount. This is evident in Fig. (□). 
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4.2 The matter power spectra 

We now discuss the effect of zt and Wf on the CDM and Q 
power spectra. A key point is that the scalar field is very light 
both before and after the transition at z t . This means that 
the associated Compton wavelength, ~ (V")^ 1 ^ 2 , of the Q 
field is very large (as in standard quintessence models) for 
all times. This implies that no clustering occurs in the Q 
field on small scales (fc > O.Ol/i Mpc -1 ). 

After the transition the potential V(Q) becomes even 
flatter (in order to obtain acceleration) and hence the Comp- 
ton wavelength increases. This means that clustering now 
only occurs on the largest scales for z < z t . This allows us 
an intuitive idea of the effect of z%. As z t is increased, we 
see that the Compton wavelength effect forces clustering to 
occur only on larger and larger scales. However, this is ac- 
tually only a fairly weak effect and the dominant variable is 

Wf. 

The effect of Wf on the Q power spectrum is straight- 
forward: for fixed zt, the closer Wf is to — 1, the less clus- 
tering occurs. Conversely, the closer Wf is to 0, the more 
clustering occurs. This is clear in the quintessence case from 
the work of Ma et al. (1999) which can be understood by 
rewriting the RHS of equation (||) in terms of the CDM 
density perturbation S c = SpcDM / PCDM which becomes 
S c [(l + Wf)pQ] 1 ^ 2 . Clearly this driving term drops to zero 
as Wf — > — 1. 

For fixed Wf, increasing zt implies that the universe 
spends less time in the dust phase where w = and hence 
the long wavelengths of SQ have less time to grow relative 
to the short wavelengths (that will not grow irrespective 
of the values of zt and Wf ). This effect is clearly visible 
in Fig. [l] where we show the global best-fit CDM power 
spectrum with z± = 1.5 together with the best-fit to just 
the LSS data which has z t = 6.5 and hence less power on 
large scales. Similarly for fixed Zt, increasing Wf towards 
zero allows more clustering on large scales relative to small 
scales. 

An important point is that simply specifying Wf = — 1 
and Zt > 1000 does not imply that the resulting model is 
the same as a ACDM model. While the transfer function 
T(k) = 1 on all scales for z < z t if Wf = — 1, this does not 
mean the initial power spectrum Pg(k) was zero, whereas 
Pa(/s) = since 8 A — by definition. This is visible in 
the C' e s of the left panel of figure (|l]). Both curves have 
fig = 0.73 and Wf = — 1. They differ due to the fluctuations 
SQ and the transition at z t = 1.5. 

This means that despite the background dynamics be- 
ing essentially equivalent for the two models for zt — > oo, 
the perturbations are not the same, and indeed one could 
consider adiabatic or isocurvature initial conditions for the 
SQ. Hence, simply showing that the recent dynamics of the 
universe favours w = — 1 does not of itself, prove that the 
acceleration comes from the cosmological constant. 

Tests sensitive to perturbations are also required. This 
may be particularly important in the case when the scalar 
field is non-minimally coupled to the spacetime curvature 
(see e.g. Perrotta & Bacciagalupi, 2001). 
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Figure 3. The redshift dependence of the luminosity distance (as 
magnitudes) minus an empty (Q = 0) universe for four different 
metamorphosis models. The z t = 0.5 model fits the data best, 
mainly due to the single data point at z = 1.7. The solid line 
effectively coincides with a ACDM model. The redshift-binned 
SNla data is from Riess et al (2001); the four dashed data points 
are experimental and were not included in the fit. 

4.3 The SNla data 

To compare with the supernovae type la measurements we 
compute the luminosity distance: 

d h {z)=c{l + z) [ du/H{u) (6) 
Jo 

from which A(m — M) — 5[logw(d,L(z)) — logio(d,L (z))] 
where dt = cz(z + 2)/(2Ho) is the empty-beam distance. 
We find that a step-function approximation to w(z) is ex- 
tremely accurate, due to the integral nature of the luminos- 
ity distance (see the dash-dotted line versus the short dashed 
line in figure (H)). 

Figure (pIT shows A(m — M) for a variety of models 
while figures (Uh and (|) show the results of the likelihood 
analysis. Models with a very recent transition from decel- 
eration to acceleration are favoured since they fit the high- 
est z supernova best, while still being consistent with the 
intermediate-^ supernovae at z ~ 0.5. Due to the large er- 
ror bars, the constraints are weak, however. Furthermore, 
since the metamorphosis models are almost indistinguish- 
able from standard ACDM models for z < z t , we have only 
constraints on zt less than the redshift of the farthest super- 
nova observed (z w 1.75). 

The dependence on the other parameters is very much 
the same as for conventional dark energy models. Not sur- 
prisingly, we recover for Qq vs wj the results of Turner and 
Riess (2001). 

4.4 BBN constraints 

Big-bang nucleosynthesis (BBN) primordial abundances de- 
pend sensitively on the expansion rate of the universe at the 
temperature T = lMeV which controls the neutron-proton 
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2D LIKELIHOODS FOR THE COMBINED DATA (CMB + LSS + SNla) 




Figure 4. The marginalised 1-d likelihood plots for our vari- 
ables (Qq , zt,w j). Left column: CMB + LSS, middle column: 
SNla, right column: Total data set. Both the CMB and SNla 
data favour a late transition (small zt) due to the corresponding 
delay in cosmic acceleration. We computed these marginalised 
likelihoods using both integration and maximisation and the re- 
sults were similar, as they should be for Gaussian likelihoods. 



ratio. Assuming that the quintessence field scales as radia- 
tion at nucleosynthesis Bean el at (2001) set the limit of 
Qq < 0.045 at 2a at T = lMeV. 

Since we assume the initial wo = 0, the scalar field scales 
as dust for z > z t and hence Qq is dynamically negligible at 
nucleosynthesis which therefore provides no constraints on 
our parameters. 

If we broadened our parameter set to include too we 
would expect BBN to set joint limits on the parameters. In 
particular, for wq = 1/3 (radiation), the BBN data would 
favour large Zt, Wf close to —1 and smaller Qq. 



5 RESULTS 

Figures (^) and (^) show our main results through the 
marginalised 1-d and 2-d likelihoods for (zt, Wf, Qq). 

We do not show the likelihoods for CMB and LSS alone 
since the LSS provides only very weak constraints on z t 
(slightly preferring higher values) and w / (no constraints at 
all). It prefers an Qq around 0.7, consistent with the CMB 
likelihood. Due to these very weak results (and the consis- 
tent result for Qq), the likelihood for the CMB data looks 
just like the one for the combined CMB+LSS data. 

The supernovae prefer (as explained in section |4~^ ) a low 
Zt, but its significance (stemming from only one supernova at 
z > 1) is too weak to change the overall likelihood by much. 
The constraints on Qq are weaker than for the other data 
sets, but consistent. Furthermore, the SNla data prefers w w 
— 1, which tightens the overall constraints on the equation 
of state somewhat, leading to w < —0.8 at the la level. 

In Figs. (|l|) and (^) we show our best fits versus the cur- 
rent CMB, LSS and SNla data and theoretical predictions 
of the standard ACDM model. 
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2D LIKELIHOODS FOR THE SNla DATA 
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2D LIKELIHOODS FOR THE CMB DATA 




0.6 0.65 0.7 0.75 0.8 



Figure 5. The marginalised 2-d likelihood plots for the com- 
bined CMB, LSS and SNla data sets showing 1 and 2<r contours 
defined as where the integral of the normalised two-dimensional 
likelihoods are equal to 0.68 and 0.95 respectively. 



The x 2 values of the overall best fit model (with Wf = 
-1.0, z t = -1.5 and Qq = 0.73) are 33 (CMB) + 36 (LSS) 
+ 4 (SNla), in total 73. On the same parameter grid, the 
best fit ACDM model has Qq = 0.73 as well. Its \ 2 values 
are 40 (CMB), 34 (LSS) and 4 (SNla), in total 78. 

We used 49 data points for the CMB, 48 for the LSS and 
7 for the supernova data. We allowed a free overall model 
normalisation plus a bias/calibration uncertainty for each 
of the three LSS and the four CMB data sets. Our phe- 
nomenological model has three free parameters, while the 
ACDM models have only one. So in total we have approx- 
imately (neglecting correlations within the experiments as 
well as between them) 49 + 48 + 7-1-4-3-3 = 93 
degrees of freedom for our model, and 95 dof for the ACDM 
models. 

We can see that both groups of models are perfectly 
consistent with current data. Given the error bars of the 
data sets, the family of ACDM models is included in our 
phenomenological models for w = — 1 and large zt- The fig- 
ures show that current data slightly prefers a low-z phase 
transition, which is still true when taking into account that 
we have to add two degrees of freedom for pure ACDM mod- 
els. On the other hand, the difference is too small to speak 
of a detection; assuming Gaussian errors and 3 "parame- 
ters of interest" (Qq, Wf, z t ), models with a A^ 2 of 5 above 
the best-fit would formally be excluded at about 83%, hence 
somewhere between 1 and 2 a. 
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6 CONCLUSIONS AND TESTS FOR LATE 
TRANSITIONS 

We have studied a phenomenological model in which the 
dark energy of the universe is described by a scalar field 
Q whose equation of state w undergoes a sudden transition 
(metamorphosis) from wo = (dust) to to/ < —0.3 at a 
specific redshift z t . 

While similar to the quintessence paradigm in prac- 
tical respects, the underlying philosophy is very different 
since we are interested in the possibility of detecting radical 
physics in the dark energy, such as the vacuum metamorpho- 
sis model (PR). We used the current CMB, large-scale struc- 
ture (LSS) and supernovae (SNla) data to constrain our 
phenomenological parameter space variables (Qq, Zt, Wf). 

The CMB and SNla data are sensitive to a transition 
if it occurs at low redshifts (z t < 3) due to the delay in the 
epoch at which cosmic acceleration can begin, relative to the 
standard ACDM models. We found that 



• The global best-fit to the current data occurs for z± = 
1.5,10/ = —1.0 and Q,q = 0.73 while the marginalised Id 
likelihood for zt peaks at z t = 2.0. The best-fit model is 
consistent with the data and is a marginally better fit than 
the best ACDM model. 

• The CMB provides the best constraints on the parame- 
ters, especially on zt due to the integrated Sachs- Wolfe effect 
(Fig. |). 



Is it possible to distinguish this metamorphosis from stan- 
dard quintessence models? 

This may be difficult since studies of quintessence favour 
similar values for iu/ and fig (Wang et al 1999, Bacciagalupi 
et al. 2002, Corasaniti & Copeland 2002) and c\ = 1 in both 
cases (Erickson et al. 2001). Traditional methods to discover 
variation in w based on the luminosity /area distance may be 
sufficient to discover transitions at z t < 2 but may be sta- 
tistically inefficient in separating metamorphosis from stan- 
dard quintessence since the distances effectively depend on 
the integral J w(z)dz (Maor et al. 2001). 

An interesting alternative (Jimenez & Loeb 2001) is 
measurements of age differences in passively-evolving galax- 
ies at different redshifts which in principle allow direct de- 
termination of w(z) and a direct test of metamorphosis if Zt 
lies in the epoch of galaxy formation (z t < 5). 

Furthermore, future very high-precision measurements 
of the CMB at 1% or better might be able to detect the per- 
turbations in the microwave background from the fluctua- 
tions in the scalar field itself, which would not be present in a 
smooth background component like a cosmological constant 
though separating this out from lensing and foreground con- 
tamination will be very difficult. 

Finally, an intriguing possibility is that the rapid tran- 
sitions studied here may provide a solution to the current 
impasse for quintessence models in explaining the varying- 
a data,wz: quintessence models can explain the apparent 
variation of a around z ~ 1 — 3 but cannot then simultane- 
ously match the results of the Okun natural reactor (Chiba 
& Khori 2001) at z ~ 0. Detailed analysis of these issues is 
left to future work. 
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